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Abstract 


Land-atmosphere feedback, by which (for example) precipitation-induced 
moisture anomalies at the land surface affect the overlying atmosphere and 
thereby the subsequent generation of precipitation, has been examined and 
quantified with many atmospheric general circulation models (AGCMs). Gen- 
erally missing from such studies, however, is an indication of the extent to 
which the simulated feedback strength is model dependent. Four modeling 
groups have recently performed a highly controlled numerical experiment 
that allows an objective inter-model comparison of land-atmosphere feedback 
strength. The experiment essentially consists of an ensemble of simulations 
in which each member simulation artificially maintains the same timeseries 
of surface prognostic variables. Differences in atmospheric behavior between 
the ensemble members then indicates the degree to which the state of the 
land surface controls atmospheric processes in that model. A comparison 
of the four sets of experimental results shows that feedback strength does 
indeed vary significantly between the AGCMs. 
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1 Introduction 


The impact of precipitation anomalies on soil moisture anomalies is self- 
evident - heavy rains induce wet soil, whereas extended dry periods induce 
dry soil. Less obvious is the impact of soil moisture anomalies on the pre- 
cipitation itself. Conceivably, a wetter soil can produce higher evaporation, 
which in turn can induce additional precipitation through both local recy- 
cling and modifications in the large-scale circulation. This land-atmosphere 
feedback, if strong, is of great interest. It may allow, for example, the trans- 
lation of soil moisture anomalies into short-term and seasonal predictions of 
precipitation. It also helps determine the climatic impacts of land-use change 
(e.g., deforestation). 

Atmospheric general circulation models (AGCMs) are popular tools for 
examining land-atmosphere feedback, largely because they include param- 
eterizations for many of the physical processes involved and because these 
parameterizations can be manipulated easily in controlled experiments. The 
list of published feedback studies is extensive [e.g., Shukla and Mintz, 1982; 
Henderson-Sellers and Gornitz , 1984; Delworth and Manabe, 1989; Oglesby 
and Erickson, 1989; Dirmeyer, 1994; and Lau and Bua, 1998; to name only 
a few]. Necessarily missing from single- AGCM experiments, however, is an 
analysis of the degree to which the experimental results are model-dependent. 
Models can certainly differ in their simulation of feedback; differences in land- 
surface parameterizations, for example, can lead to differences in the response 
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of evaporation to precipitation anomalies, and differences in boundary layer 
and convection parameterizations can lead to differences in the atmosphere’s 
response to surface evaporation and sensible heat flux. Clearly, in any feed- 
back study, an evaluation of simulated feedback strength against observations 
is desirable. 

Unfortunately, while a few indirect approaches have been employed at 
the regional scale [e.g. , Findell and Eltahir, 1997], the direct quantification 
of real-world feedback strength at the global scale from available observations 
is extremely difficult, if not impossible. The validation of simulated feedback 
strength is indeed beyond the scope of this paper. This paper instead focuses 
on a lesser, but still very important, aspect of the problem: the extent to 
which simulated feedback strength varies between different AGCMs. This 
variation is a measure of the degree to which feedback-related model results 
are model dependent. In a sense, it measures the uncertainty inherent in our 
understanding of feedback and our ability to model it. 

The problem of intermodel variation in feedback strength is addressed 
here by having several AGCMs perform the same, highly controlled nu- 
merical experiment. Four models participated in the intercomparison: the 
NASA Seasonal-to-Interannual Prediction Project (NSIPP) AGCM, the Cen- 
ter for Ocean-Land-Atmosphere Studies (COLA) AGCM, the National Cen- 
ter for Atmospheric Research Community Climate Model Version 3 cou- 
pled to the Biosphere- Atmosphere Transfer Scheme (CCM3/BATS), and the 


3 



United Kingdom Meteorological Office AGCM (HadAM3). The experiment 
is described in section 2, and the results of the intercomparison are presented 
in section 3. Section 4 offers some interpretation of the results. 

2 Experiment Design 

The design of the experiment is illustrated in Figure 1. The experiment has 
two parts. In the first part, the AGCM, fully coupled to its own land surface 
model (LSM) but forced by prescribed sea surface temperatures (SSTs), is 
run over a selected month. At each time step in this simulation (hereafter 
labeled simulation W, for “write”), the values of all land surface prognos- 
tic variables at every grid cell are recorded into a special data file. The 
recorded prognostic variables include soil moisture contents at all vertical 
levels, temperatures at all vertical levels, canopy interception reservoir con- 
tent, and various variables characterizing snow, if snow is present. The one- 
month experiment is then repeated 15 more times, using 15 different sets 
of atmospheric and land surface initial conditions, to obtain an ensemble of 
16 one-month simulations (simulations W2-W16). The prognostic variables, 
however, are only recorded during Simulation Wl. 

The second part of the experiment consists of another 16-member ensem- 
ble of one-month simulations, using solar forcing for the same month as before 
and using the same prescribed SSTs. Again, the ensemble members use dif- 
ferent atmospheric initial conditions. At every time step of every simulation, 
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the updated values of all land surface prognostic variables are discarded and 
then replaced by the corresponding values for that time step from the data 
file written in Simulation Wl. Thus, in this ensemble, all member simula- 
tions (simulations R1-R16, where R denotes “read”) are forced to maintain 
precisely the same time series of (geographically-varying) land surface states. 

We focus mainly on precipitation in this study. The idea is simple: we 
examine the degree to which the time series of precipitation rates in simu- 
lations R1-R16 are similar. If they are similar, even after “subtracting out” 
the effects of SSTs and other intramonthly signals through an analysis of 
simulations W1-W16, then we can say that the evolution of the atmosphere 
is indeed largely governed by land surface conditions. The degree to which 
the intra-ensemble similarity differs between GCMs provides one measure of 
how land-atmosphere feedback strength varies between models. 

The set-up of the experiment did differ slightly between the four partic- 
ipating GCMs. For example, for various reasons, each GCM used its own 
set of prescribed SSTs, though each was derived from datasets used in vari- 
ous phases of the Atmospheric Model Intercomparison Project [Gates, 1992]. 
These differences, which should not affect the basic outcome of the experi- 
ment, are outlined in Table 1. 
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3 Results 


To quantify the degree of “similarity” in the time series of precipitation 
amongst the members of an ensemble, we employ an approach used by Koster 
et al. [2000]. First, we choose an aggregation period. In the analysis below, 
we look at 3-day totals of precipitation, P; each July simulation thus provides 
a times series of ten 3-day totals. We then compute an ensemble average 
time-series, P. For each time period n, we compute 

A -is E** (1 > 

iU t=i 

where i loops over the 16 ensemble members. Next, we compute the variance, 
0-2 j 0 f p across all ensemble members and time periods (i.e., across 160 
values) and the variance, <7p, of P across all time periods (i.e., across 10 
values). This allows us to compute Op, a measure of time series similarity: 



Note that if each ensemble member produces exactly the same time series of 
P, then dp will equal dp, and ftp will equal 1. If, however, the time series are 
completely uncorrelated, then d 2 p will be approximately tfp/16, and 0 P will 
be about 0. Thus, 0 P varies from 0 to 1, with values closer to 1 indicating a 

greater degree of precipitation similarity. 

Op essentially measures the ratio of the signal variance to the total vari- 
ance. (A similar diagnostic was suggested by Rowell et al. [1995]). Figure 
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2 illustrates the nature of ftp graphically. The top plot in the figure shows 
the time series of precipitation at a specific grid cell for each of the 16 sim- 
ulations in the NSIPP R ensemble. Note that the precipitation is low for 
the month until day 20, when it shoots up for each simulation. The obvious 
coherence between the different time series is reflected in the high ftp value 
(0.85) at this grid cell. In the bottom plot, which shows the 16 time series 
at a different grid cell, the coherence is absent - the precipitation generated 
in one simulation is essentially independent of that in any other simulation. 
Precipitation at this cell is thus controlled by chaotic atmospheric dynamics 
rather than by SST or land conditions, ftp for this grid cell is very low (0.07). 

To relate ftp to land control, we must properly account for seasonal vari- 
ations in SSTs and anything else outside of the land surface that can induce 
intra-monthly trends in precipitation. We do this simply by calculating ftp 
separately for the R and W ensembles. flp(R) represents the similarity in 
precipitation induced by all factors, including the specified set of land states. 
ftp(VF) represents the similarity induced by everything but the specified land 
states. The difference ftp(R) — ftp (VP) thus gives a first-order indication of 
the land’s impact on the evolution of the atmosphere. 

Global maps of fl P (R) — fl P (W) are provided in Figure 3 for all four 
AGCMs. The salient result is a wide disparity in the diagnostic between the 
models. The impact of land conditions on atmospheric processes is clearly 
largest for the NSIPP model. The COLA and CCM3/BATS models have 
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similar Op(f?) — Op (kb) distributions, with values of 0.2 or less almost ev- 
erywhere, and the HadAM3 model has what appears to be an even weaker 
land-atmosphere connection. 

Because the choice of a 3-day averaging period was somewhat arbitrary, 
other averaging periods were examined as well. The relative behavior of the 
models is similar when the precipitation is averaged over 1-day and 6-day 
periods (not shown). In general, however, Q P (R) — Q P (W) for a given model 
increases as the precipitation aggregation period increases. - 

4 Discussion 

4.1 Coherence of Surface Fluxes 

One potential explanation for low values of Op(f?) — Op(kb) in Figure 3 in- 
volves the response of the surface turbulent fluxes to the imposed surface 
states. Due to variations in the atmospheric forcing, the time series of evap- 
oration or sensible heat flux amongst the members of an ensemble may not 
look the same, even given identical time series of soil moisture and temper- 
ature. If the time series of turbulent fluxes were not the same, and if the 
effect of the land surface on precipitation were mostly through these fluxes, 
then the derived values of Q P (R) - Op (IF) would necessarily be low. 

A warning is in order here — the response of precipitation to land surface 
conditions may be more complex than this. For example, the temperature 
and humidity conditions in the atmospheric boundary layer, which help de- 
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termine precipitation, may be strongly guided by surface temperature and 
moisture states. Although the surface fluxes are responsible for communi- 
cating the surface states to the boundary layer, the state of the boundary 
layer may nevertheless correlate better with the surface states themselves 
than with the time-integrated surface fluxes examined in this section. 

With this caveat, we analyze the time series of evaporation rates produced 
by the GCMs by defining, in direct analogy to f2p, the diagnostic Q#: 


n 


E 


16<4 - a% 
1 


( 3 ) 


where <r| is the variance of evaporation across all ensemble members and cr| 
is the variance of evaporation in the ensemble mean time series. £Ie thus 
measures the degree of similarity among the time series of evaporation rates 
(and loosely, via energy balance considerations, among the time series of 
sensible heat fluxes) produced by the different members of an ensemble. As 
with Clp, a value of 1 implies that all 16 time series are identical, whereas a 
value of zero implies that the 16 time series are completely uncorrelated. 

Figure 4 shows the global distribution of Qe{R) ~ for each of the 

AGCMs. In each case, Qp is high over much of the globe. Indeed, in many 
regions, specifying surface moisture and temperature states in this experi- 
ment is roughly equivalent to specifying time series of the surface turbulent 


fluxes. If fip is correspondingly low in these regions, as it is in particular for 
the COLA, CCM3/BATS, or HadAM3 GCMs, then we can conclude that 
the modeled atmosphere in these regions does not respond strongly to the 
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local surface fluxes. That is, a low ftp value in the presence of a high Op 
value is probably best explained by the model’s atmospheric formulations, 
presumably those of boundary layer processes and moist convection, though 
it might also (or instead) result from a low Qe value in a critical remote 
region. 

In some regions, flp is indeed fairly low. Both the CCM3/BATS and 
the HadAM3 models, for example, show low values over a large fraction 
of northeastern Asia, and the HadAM3 model also show r s low values over 
much of North America. Such findings for HadAM3 are consistent with 
Gedney et aids [2000] suggestion that evaporation in the HadAM3 model 
is rarely moisture limited. If the land surface’s effect on precipitation is 
mostly through the surface fluxes, and if this effect is mostly local, then 
Figure 4 does identify some regions for which the land-atmosphere coupling 
is weakened by an insensitivity of surface moisture and energy fluxes to the 
surface prognostic states. 

4.2 Idealized Nature of Experiment 

One curious aspect of Figure 3 is the geographical structure of the NSIPP 
model’s f l P distribution, which differs somewhat from the distribution of sea- 
sonal f Ip values derived by Koster et al. [2000] with the same AGCM. Koster 
et al. [2000] place the largest Q P values in the transition zones between arid 
and humid regions, wdiereas Figure 3 places them in somewhat different ar- 
eas, farther to the north, for example, in North America. The difference 
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presumably results from the different designs of the experiments. Koster et 
al. [2000] effectively prescribe only land surface soil moisture contents. This 
has direct relevance to predictability studies, since soil moistures may be 
predictable (and thus effectively “specified”) on seasonal timescales due to 
their inherent memory. The present experiment prescribes diurnal variations 
of both surface moisture conditions and surface temperatures. Apparently 
the added specification of temperature leads to a significant increase in Qp 
in certain areas. 

This last point underscores the idealized nature of the present experiment. 
The experiment does not address the long-term predictability of precipitation 
in the participating models. It simply allows a highly controlled comparison 
of the models’ coupling behavior - a comparison of the degree to which the 
generation of precipitation is guided by the full land surface state rather 
than by chaotic atmospheric dynamics or external forcings (SSTs). Clearly 
the models disagree about the strength of land-atmosphere coupling. This 
uncertainty is often not acknowledged in land-atmosphere feedback studies. 

We must emphasize again that the proper level of land-atmosphere feed- 
back — the proper range of fip (or even 0^) — is simply not known. This 
analysis makes no attempt to state which model’s coupling behavior is the 
most realistic. For a proper evaluation, long-term observational studies of 
boundary layer behaviof are required. Furthermore, the low f Ip values in 
Figure 3 for the COLA, CCM3/BATS, and HadAM3 models do not imply 
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that these models have a low potential for predictability. The COLA group 
recently performed a series of experiments similar to those of Koster et al. 
[2000]. An analysis of these experiments shows that the coupling behavior 
of the COLA model on seasonal timescales is similar to that of the NSIPP 
model — in effect, flp for both models is high over a significant fraction of 
the earth. Apparently the noise behind the COLA model’s low Qp values in 
Figure 3 is smoothed out at much longer timescales. 

5 Summary 

A highly-idealized, simple, and inexpensive AGCM experiment has been de- 
vised that illustrates some key aspects of simulated land-atmosphere coupling 
behavior. Four AGCM groups have performed the experiment for the present 
study, and a comparison of their results shows that the apparent strength 
of land-atmosphere coupling on short (daily to weekly) timescales does vary 
significantly between the models. A strict evaluation against observations 
is not provided, since the strength of coupling in the real world is difficult, 
if not impossible, to measure. Nevertheless, the intermodel differences are 
important in themselves, since they illustrate the uncertainty with which we 
understand the various processes (particularly atmospheric boundary layer 
and convection processes) that control the coupling. Further studies of cou- 
pling strength are needed, particularly with field observations, given the cen- 
tral role played by land-atmopshere feedback in many published and ongoing 
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numerical studies of climate variability and predictability. 
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Figure Captions 

Fig. 1 Basic design of the experiment, as performed by all participating mod- 
els. 

Fig. 2 Superposed time series of precipitation for two grid cells: one in which 
rip is high (top), and one in which Q P is low (bottom). 

Fig. 3 Global fields of Qp(i?) — fip(fF), as generated by each of the partici- 
pating AGCMs. 

Fig. 4 Global fields of ftp(i?) — rip (IT), as generated by each of the partici- 
pating AGCMs. The values of this diagnostic over the ocean were not 
available for HadAM3. 
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Model 

SSTs Used 

Surface 
Grid Res. 

Key References for AGCM 

Notes 

NSIPP 

July 1988 

4° (long.) 
x 5° (lat) 

Dynamics: Suarez and 
Takacs [1995] 

Land: Koster and 
Suarez [1992, 1996] 
Convection: Moorthi and 
Suarez [1992] 

Initial atmosphere 
and land states for 
W and R from July 1 
1988 conditions in 16 
parallel AMIP runs. 

COLA 

July 1986 

1.8°x 2.8° 

Physics: Kinter et al. [1997] 
Dynamics: Kiehl et al. [1998] 
Land: Xue et al. [1991], 
Dirmeyer et al. [1999] 
Convection: DeWitt 
et al. [1996] 

Full description of 
experiments in 
Dirmeyer et al. 
[2001], 

CCM3/ 

BATS 

July 1983 

j 

i 

2.8° x 2.8° 

GCM: Kiehl et al. [1998] 
Land: Dickinson et al. [1993] 
GCM-Land comb.: Hahmann 
and Dickinson [2001] 

W2-W16 began June 
15 (first 15 days 
thrown out); June 15 
land states were 
identical; June 15 
atmosphere states 
from 16 sequential 
days. 

HadAM3 

July 1981 

3.75°x 2.5° 

General: Pope et al. [2000] 
Land: Cox et al. [1999] 
Convection: Gregory and 
Rowntree [1990] 

Initial atmosphere 
and land states for 
W from a previous 
AMIP run: July 1 
1979-1994 inclusive. 
Prescribed land states 
for R from 1981 run, 
with land states 
initialized to 
climatology. 


Table 1: Models used in the experiment and details of the experimental 
set-up. 
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PART 1: ESTABLISH A TIME SERIES OF SURFACE CONDITIONS (Simulation W1). 



(Repeat without writing to obtain simulations W2-W16.) 


PART 2: RUN 16-MEMBER ENSEMBLE, WITH EACH MEMBER FORCED TO 
MAINTAIN THE SAME TIME SERIES OF SURFACE PROGNOSTIC VARIABLES 
(Simulations R1-R16). 



Figure 1: Basic design of the experiment, as performed by all participating 
models. 
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Figure 2: Superposed time series of precipitation for two grid cells: one in 
which f2p is high (top), and one in winch Qp is low (bottom). 




Figure 3: Global fields of Q P (R) - Q P (W), as generated by each of the 
participating AGCMs. ^ 






Figure 4: Global fields of Q E {R) ~ ^(PF), as generated by each of the 
participating AGCMs. ^ 








Popular Summary: 


'‘Comparing the Degree of Land-Atmosphere Interaction in Four Atmospheric Gen- 
eral Circulation Models 7 ’, by R. Koster. P. Dirmeyer, A. Hahmann, R. Ijpelaar, L. 
Tyahla, P. Cox, and M. Suarez 

A wet soil that results from a rainstorm may lead to higher evaporation rates 
in the weeks to follow, and this higher evaporation may induce additional precipi- 
tation, either through local moisture recycling or through changes in the large-scale 
atmospheric circulation. Similarly, a drought may lead to dry soil and thus low 
evaporation rates, and the low evaporation may help stifle subsequent precipitation, 
thereby prolonging the drought. These are examples of land-atmosphere feedback, 
by which atmosphere-induced land surface anomalies in turn affect the atmosphere. 
Land-atmosphere feedback is extremely difficult, and often impossible, to quantify in 
the real world. Given its importance, though, it has been examined extensively with 
many atmospheric general circulation models (AGCMs). 

Unfortunately, many of these studies give conflicting results. This leads to some 
obvious questions: (1) To what extent is the strength of simulated feedback model 
dependent? (2) How can we quantify intermodel differences in feedback strength in 
an objective way? In this paper, we address these questions with a carefully-devised 
AGGM experiment, an experiment with two very important characteristics. First, 
it directly quantifies land- atmosphere feedback in the AGCM. This is achieved by 
imposing the same land surface states in each member of a 16-member ensemble and 
then examining the similarity of the precipitation generated by the different ensemble 
members. Second, the experiment is simple enough and computationally inexpensive 
enough to be accessible to many different modeling groups. Indeed, for the present 
paper, four groups performed the experiment. 

The models do show substantial differences in the simulated strength of land- 
atmosphere feedback. Such differences, objectively determined here for the first time, 
underline the uncertainty inherent in our understanding of land impacts on climate. 
The four-model comparison demonstrates the viability of the intercomparison exper- 
iment. which we hope will expand to include a great many more models in the near 
future. 
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